Nonlocality in mesoscopic Josephson junctions with strip geometry 
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We study the current in a clean superconductor-normal-metal-superconductor junction of length 
d and width w in the presence of an applied magnetic field H. We show that both the geometrical 
pattern of the current density and the critical current I c (&) , where $ is the total flux in the junction, 
depend on the ratio of the Josephson vortex distance a — &o/Hd and the range r ~ y/d£ N of the 
nonlocal electrodynamics ($o = hc/2e, £ N = Tiv F /2iiT, and r(T — » 0) ~ d). In particular, the 
critical current has the periodicity of the superconducting flux quantum $o only for r < a and 
becomes, due to boundary effects, 2<E>o (pseudo-) periodic for strong nonlocality, r > a . Comparing 
our results to recent experiments of Heida et al. [Phys. Rev. B 57, R5618 (1998)] we find good 
agreement. 



PACS numbers: 74.50. +r, 74.80.Fp 

Quantum mechanical interference effects render the 
electrodynamics of mesoscopic samples nonlocal. In par- 
ticular, nonlocality is a key element entering the under- 
standing of the magnetic response and the transport in 
SNS-junctions (s-wave superconductor-normal-metal-s- 
wave superconductor junctions) and SN-proximity sand- 
wiches. The strength and relevance of the nonlocality 
in general depends on the dimensions of the system, the 
normal- metal coherence length £ w = Tiv F /2irT, and the 
elastic scattering length I p]-p|. In this paper we show, 
how the different length scales enter the magnetotrans- 
port problem of a mesoscopic SNS-junction to produce a 
shift in the (pseudo-) periodicity of the critical current 
from $0 to 2$ . 

After the discovery of the Josephson effect in SIS 
tunnel junctions 0, interest turned to metallic links 
of the SNS type |g], where the current is conveniently 
described in terms of Andreev states trapped within 
the normal-metal region |fj. For a wide junction, the 
current density and the supercurrent in the presence 
of a magnetic field H have been calculated by Antsy- 
gina et al. 0, who found a $o periodicity in the crit- 
ical current. Continuous progress in nanofabrication 
technology made it possible to investigate mesoscopic 
superconductor-semiconductor heterostructures, see Ref. 
H for a study of the fluctuations in the critical cur- 
rent and its quantization in a superconducting quan- 
tum point contact. Recently, Heida et al. H, investigat- 
ing S-2DEG-S-junctions (s-wave superconductor-two- 
dimcnsional-electron-gas-s-wave superconductor junc- 
tions) of comparable width w and length d, have mea- 
sured a 24>o periodicity of the critical current instead of 
the standard $o- A first attempt to explain this find- 
ing is due to Barzykin and Zagoskin fllO| ] ; considering the 
point-contact geometry of Fig. 1(a) with open boundary 
condition in the metal, they indeed recover a 2<I>o peri- 
odicity in the limit w/d — > 0. However, the experiment 
of Heida et al. M is carried out in the strip geometry 



of Fig. 1(b) and involves dimensions w ~ d of the same 
order. In the present paper, we determine the critical 
current I c through a clean SNS-junction in the presence 
of an applied magnetic field H, taking proper account of 
the reflecting boundaries in the normal-metal character- 
istic for the strip geometry of Fig. |l|( 
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FIG. 1. (a) Junction with a point-contact (open) geom- 
etry as discussed in [10], where w is the width of the two 
superconductors, (b) The junction studied here has a strip 
geometry with w the width of the normal conductor, (c) The 
magnetic field H is applied in the z-direction and the coordi- 
nate system is chosen symmetric with respect to the junction 
center. The current density in the point P involves contri- 
butions from all trajectories P parameterized by the angle if. 
(d) The phase difference 7 along the trajectory F can be ex- 
pressed through the enclosed flux <f>; using the trajectory To 
as our reference, the flux through the areas above (below) Po 
contributes with a positive (negative) sign. 



We find that the periodicity of the critical current 
changes from $0 to 2$o as the flux through the junc- 
tion increases. At low temperatures the crossover to the 
2<&o periodic current appears at a flux ~ &Qw/d, thus 
explaining the result of Heida et al. M, who found a 2$ 



periodic pattern for all fields in devices with w/d ~ 1. 

In our derivation, we neglect screening effects by the 
induced supercurrent, which is justified as long as H > 
y/8&ojc/cd', where j c denotes the critical current density 
of the junction and d' — d + 2A with A the penetration 
depth of the two bulk superconductors fllj| . Since all 
the length scales in our system (the dimensions w and d, 
the normal- metal coherence length £ N ) are much larger 
than the Fermi wave length A^, we base our calculations 
on the Eilenberger equations [Q for the quasi-classical 
Green functions and extract the current density in the 
standard way M. 

The SNS-junction we study is sketched in Fig. [y. In the 
quasi-classical formulation, the current density in a point 
P results from contributions over all trajectories of quasi- 
particles connecting one interface to the other through 
P. In a junction of infinite width, all trajectories go 
straight through the junction. In the case of a finite junc- 
tion, boundary conditions at the normal-metal-vacuum 
boundary have to be applied, which we idealize through 
the assumption of specular reflections. Furthermore, we 
adopt the usual approximations: perfect Andreev reflec- 
tions at the SN-interfaces and a coherence length £o in the 
two superconductors with £o -C d, allowing for a step-like 
approximation of the order parameter A RMO]. The 
quasi-classical Green function is calculated by matching 
the partial solutions in N and S at the interfaces. For 
the current density j, we arrive at a generalization of the 
results of Antsygina et al. J7J . Explicitely, for finite tem- 
peratures with d ^> £ N , j takes the form, 
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where p = (cos(y),sin(<^),0), l(ip) — d/ cos(ip) is the 
length of a trajectory with slope tp, and the critical cur- 
rent densities are 



3c,t = PJc,o exp 



& 



and j c 



ne 1 $o 
mc 2d 



(3) 



In (0) , n denotes the electron density in the normal con- 
ductor and p sw 12/vr for T < T c , lim T _ Tc p~l- T/T c . 
While in the low temperature limit all harmonics sin(k"f) 
(k = 1,2, . . .) contribute to the current density |13|, at 
finite temperatures, thermal smearing of the Andreev 
levels leads to a suppression of the higher harmonics 
oc exp(— kd/^ N ) and only the first term oc sin(7) sur- 
vives. An individual trajectory contributes with a weight 
oc exp(— Z/£jv) at finite- and oc d/l in the low temperature 
limit. In a wide junction, 7 takes the form [J7|, 



2tt 

j(x,y;ip)=jo-—Hd'{y-xta,n(ip)}. (4) 

The more general expression derived here results in the 
gauge invariant phase difference 



2ir f 
7(0:, y; <p) = Aip - — / A • ds, 



(5) 



where Aip denotes the phase difference between the two 
superconductors and T is the path which goes through 
the point (x, y) with slope ip. Combining the current ex- 
pressions (Q) or (|2|) and (g) with the Maxwell equation 
V 2 A = — 47rj[A, Aip]/c, we obtain the transverse vector- 
potential A and thus can solve the full screening problem; 
in the case of a tunnel junction (where the trajectories 
are reduced to the one with ip = 0), numeric and analytic 
calculations have been given by Owen and Scalapino JL4| . 
In the following, we neglect screening and concentrate 
on junctions with the strip geometry of Fig. |lj(b), in- 
cluding the (reflecting) trajectories T in (|l|) and (g). Wc 
express the gauge invariant phase difference (J5|) in terms 
of the flux <fi enclosed by T and the reference path To and 
obtain, 



i{x,y,<p) = 70- 



2ir(f>{x 1 y;ip) 



(6) 



where for negligible screening 4>(x,y;ip) = HA(x 1 y;p) 
and A is the properly weighted enclosed area, see 
Fig. |l](d) . The area A is calculated as a function of the 
number of reflections the trajectory T undergoes (in the 
following called the 'order' of the trajectory). The point- 
contact geometry of Fig. |l](a) then is described by the 
order-zero trajectories alone |IG], while in the strip ge- 
ometry of Fig. pl(b), higher orders have to be included. 
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FIG. 2. The 'bow ties' of width r display the ensemble 
of trajectories contributing to the current density in a given 
point. The arrows indicate the current flow, (a) For weak 
nonlocality, r < a , the current flows straight through the 
junction, (b) The strong nonlocality, r > a , leads to the 
formation of vortex-like domains of circular flow. 

The geometrical pattern in the current density j de- 
pends strongly on the sample dimensions d and w, the 
normal- metal coherence length £ N , and the applied field 
H . At finite temperature, the current density in P draws 
its weight from trajectories with ip < y/£ N /d, allowing us 
to introduce the transverse nonlocality range r = \/£, N d 



(in the low temperature limit, <p ~ 1 and we define r = d). 
This range of nonlocality has to be compared to the scale 
a a = $o/Hd' of transverse variations in j (see Fig. g): For 
weak nonlocality, r < a , the flow is uniform along x with 
amplitude j c and changes direction on a distance a /2 go- 
ing up the y-axis. This contrasts with the strongly nonlo- 
cal case r > a„, where the current density forms domains 
of left and right going circular flow. While the local case 
is similar to that in a tunnel junction, the pattern in the 
norilocal situation reminds of the usual vortex structure 
in a superconductor, see Fig. 0. Explicitely, for finite 
temperatures with ii > („, the current density of the 
order-zero trajectories is given by 



3x(x,y) 

3c,T 

3y(x,y) 

3c,T 

where 



- sin ( 70 j exp [-a(*)] , 

2ny\ 6v 2ttx 

cos 70 — exp -ai , (7) 

a I d a 



a(x) 




d\ (2-kxY 



(8) 



For weak nonlocality, a < r, the exponent remains ap- 
proximately constant in the normal part, a(x) ~ ct(0), 
leading to a uniform current flow, while for strong nonlo- 
cality, a > r, a(x) grows as x approaches the interfaces, 
a(±d/2) 3> a(0), such that the current concentrates in 
the middle of the junction. For a > r, the higher order 
trajectories lead to a refinement of the current pattern, 
see Fig. || (similar results are obtained in the low tem- 
perature limit). 



versus flux <f> = Hd'w in the junction. In the case of 
weak nonlocality, r < a , the relevant contribution to 
the critical current comes from the order-zero trajecto- 
ries resulting in a $o periodicity. For strong nonlocality, 
r > a , higher orders are relevant and lift the order-zero 
result as shown in Fig. — the periodicity of the criti- 
cal current changes to 2$o- To be specific, we discuss in 
detail the orders 0, 1, and 2 for the case of finite tempera- 
tures with d 3> £ N (the qualitative arguments for d <C £, N 
are similar). 
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FIG. 4. The critical current for d/£ N — 5 and w/d = 1/3. 
The solid curve shows the full critical current and the dashed 
curves are the contributions from the orders 0, 1, and 2. The 
orders and 1 oscillate with periodicity 2<I>o, while the sec- 
ond order decreases monotonically, remaining always positive. 
The current pattern produced by the orders and 1 is lifted 
by the order 2 contributions, and the critical current attains 
the periodicity 2$o- 
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FIG. 3. The current density for d/£, N = 3, r/a « 3.5, and 
$/$o = 3. Apart from the vortex-like domains in the middle 
produced by order-zero trajectories, additional circular flow 
is set up by the order-one trajectories. 

The ratio r/a and its associated characteristic current 
pattern manifest themselves in the (pseudo-) periodicity 
of the critical current, 



i c m 



max 

70 



uj/2 



w/2 



dy j x (0,y;7o,$), 



(9) 



For w > r > a , the critical current due to the order- 



zero trajectories takes the form 15 



J(°>(*) 



V2I c ,t w cos(7r<i>/$o) 
0F r (7T$/$ ) 2 



(10) 



where I c j> = wj C T- F° r the first-order trajectories, we 
numerically find again a 2<& (pseudo-) periodic contri- 
bution. Both components vanish with field oc l/$ 2 . 
The second- and all following even-order trajectories ex- 
hibit a large current amplitude of order j c on a scale 
a„ oc 1/$ in the junction center (0,0), a consequence of 
the (^-independence of the gauge invariant phase differ- 
ence 7 along trajectories through (0,0). Their contri- 
bution scales oc 1/$ and therefore dominates over the 
zcroth- and first-order terms at large enough fields — as 
the strongly nonlocal limit with a < r is reached, the 
periodicity changes over to 2<& . Samples with a small 
width w < r are always in the strongly nonlocal limit and 
their current pattern is 2$o periodic throughout the en- 
tire field axis. At low temperatures, the condition w < r 
transforms into the geometric requirement w < d. In 



the very limit w/d — > the periodic modulation disap- 
pears and the solution goes over into the zero field result, 
I c (&) —> lc(0)- The complete classification is given in Ta- 
ble I. 




FIG. 5. The critical current in the low temperature limit, 
d <C £jv, for w/d — 0.8, 0.9, 1.1, 1.5 (from top to bottom). 
Successive plots have been offset by 0.35. Measured data ([9], 
diamonds) is shown for the case w/d = 0.9. 

Recently, Heida et al. |J observed such a 2$o peri- 
odicity in strip like (w ~ d) S-2DEG-S junctions made 
from Nb electrodes in contact with InAs operating at low 
temperatures T = 0.1 K (similar junctions have been con- 
structed by Takayanagi et al, see Ref. [||). As the total 
flux through the junction is difficult to determine in the 
experiment, Heida et al. had to infer their 2$o periodic 
structure from a fit on four samples with different ra- 
tios w/d ranging from 0.9 to 2.2. In Fig. [|, we present 
the results of our numeric calculations for the strip ge- 
ometry, where we have properly taken into account the 
finite penetration depth of the flux into the supercon- 
ducting banks. While geometries with w/d < 1 clearly 
exhibit a 2$ periodicity throughout the entire field re- 
gion, a $ -component starts to develop at low fields in 
wide junctions. The comparison with the data of Heida 
et al. (w/d — 0.9) gives a satisfactory description of the 
pseudo-periodic structure. 

In conclusion, we have demonstrated that the current 
density and the critical current in a clean SNS-junction 
with strip geometry depend crucially on the ratio r/a 
between the nonlocality range r and the vortex distance 
a . The period of the critical current depends not only 
on the dimensions of the junction and the normal-metal 
coherence length £ N , as it is the case in a point-contact 
geometry ]!(]], but also on the applied magnetic field H. 
In particular, we obtain a 2<E>o periodicity in the whole 
current pattern at w ~ d, whereas for a point contact 
geometry, the 2$ periodicity is reached only in the limit 



w/d — ► [Q . The numerical results for the strip geom- 
etry are in agreement with the experiment of Heida et 
al. ||. For wider junctions, we predict a crossover from 
a $o- to a 2$o-periodicity at high fields. 

We thank D. Agterberg and V. Geshkcnbein for stim- 
ulating discussions throughout this work. 
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w/d — > 0, where the order-zero critical current is given 

by the result conjectured in Ref. M. 



TABLE I. The periodicity of the critical current is con- 
trolled by the three parameters w/d, w/r, and r/a . The 
table has to be read as a flow chart, starting at the top row 
and selecting the proper condition proceeding down the rows. 
Note that always w/d < w/r. 



ratio 


value 


w/d 
w/r 
r/a 


> 1 

> 
< 1 


< 

1 
> 1 


1 
< 1 


-> 


period 


$o 


2$o 


none 



